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We study the problem of mapping an unknown mixed quantum state onto a known pure state 
without the use of unitary transformations. This is achieved with the help of sequential measure- 
ments of two non-commuting observables only. We show that the overall success probability is 
maximized in the case of measuring two observables whose eigenstates define mutually unbiased 
bases. We find that for this optimal case the success probability quickly converges to unity as the 
number of measurement processes increases and that it is almost independent of the initial state. In 
particular, we show that to guarantee a success probability close to one the number of consecutive 
measurements must be larger than the dimension of the Hilbert space. We connect these results to 
quantum copying, quantum deleting and entanglement generation. 



I. INTRODUCTION 



During the last two decades a major research effort has been conducted in the emerging field of quantum information 
theory jlj . Much of this activity starts with the observation that the capacity of physical systems to process, store 
and transmit information depends on thir classical or quantum nature [2|. Quantum algorithms, that is, algorithms 
based on the laws of Quantum Mechanics, show an enhancement of information processing capabilities over their 
classical counterparts. A large collection of quantum communication protocols such as quantum teleportation 0], 
entanglement swapping Q , quantum cloning 0, [|| and quantum erasing reveal new forms of transmitting and 
storing classical and quantum information. Most of these protocols have already been experimentally implemented 
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A common assumption concerning quantum algorithms and quantum communication protocols is the capacity of 
performing transformations belonging to a fixed but arbitrary set of unitary transformations together with measure- 
ments on a given basis. In this article we study the problem of mapping a mixed initial state onto a known pure state 
using measurements as the only allowed resource, that is, a measurement driven quantum evolution. We show how 
this problem connects naturally to generation of quantum copies, quantum deleting and entangled states generation. 

This article is organized as follows: in section [H] we study the problem considering states belonging to a two- 
dimensional Hilbert space. In section ITTT1 we generalize to the case of a d-dimcnsional Hilbert space and show that 
mutually unbiased bases optimize the overall success probability. Section llVI presents the case of m target states in a 
d-dimensional Hilbert space. In section[V]we summarize our results. 



II. TWO-DIMENSIONAL CASE 



Let us consider a quantum system described by a two-dimensional Hilbert space. Initially, the system is in a mixed 
state p. Our goal consists in mapping this state onto the known target state \ip) by using quantum measurements as 
the only allowed resource. 

In order to accomplish this task we define a non-degenerate observable (p. Its spectral decomposition is 

= %)<p|+Aj>_l)(p_l| > (1) 

where the \tp) and \<p±) states are eigenstates of (p with eigenvalues A and Aj_ respectively. Thereby, the target state 
must belong to the spectral decomposition. 

A measurement of the (p observable onto the p state projects the system to the target state \(p) with probability 
p = ((p\p\ip). In this case we succeed and no further action is required. However, the process fails with probability 
1 — p when the measurement projects the system onto the \(p±) state. Since this state cannot be projected to \tp), 
the target state, by means of another measurement of if, it is necessary to introduce a second observable 9 whose 
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nondegenerate eigenstates |0) and |1) are 

|0) = cos(0)|^) W*sin(0)|^x) 

11} = -e-^ S in(0)|^>+cos(0)|^), (2) 

with 9 and <f> being real numbers. 

A measurement of 9 projects the \p>±) state onto the state |0) or Since both states have a component on the \p>) 
state, a second measurement of (p allows us to project again, with a certain probability, to the target state \<p). The 
probability p' that this procedure fails after a first measurement of tp but is successful after a consecutive measurement 
of the 9 and <p operators is 

p> = (^|p|^)(|(0|^)| 2 |M0)| 2 + |<1>j.)| 2 |(^|1)| 2 ) = \{<p x \p\<P±) s™ 2 {29), (3) 

where we have used Eq.Q. Then, the success probability in the sequence of measurements (M (ip) M (9)) M (ip) is 

p + p' = 1- (vjj.IpIvj.X1 - isin 2 (20)). (4) 

Similarly, the success probability p s of mapping the initial state p onto \<p), the target state, after applying the 
consecutive measurement proceses [M (ip) M (9)} N M (ip) , that is, a measurement of ip followed by N measurement 
processes each one composed of 9 followed by <p, is given by 

p StN = i-(v ± \ P \ v± )li-J2\(j\v±)\ 2 Mfi\ 2 ] > (5) 

V 3=0 



or equivalently 



i>,..\ I Ivj l/'lv±) ( l-^sin 2 (20)) . (6) 



Clearly, the extreme values 9 — and 9 = tt/2 correspond to observables ip and 9 defining the same basis. Conse- 
quently, in this case the success probability becomes simply (tp\p\(p). The expression © indicates that the success 
probability p Sl N can be maximized by choosing 9 — n/4. In this case we obtain 

_ 1 _ (yjjplyx) (7 , 
Fig. JIJi) illustrates the behavior of the maximum success probability Ps, max. 

as a function of N for different values of 
((p\p\ip). We observe that p s ,max quickly converges to 1 almost independently of the (ip\p\(p) even if the initial state p 
belongs to a subspace orthogonal to \<p). For instance, in this particular case, if the success probability for the first 
measurement of (p vanishes, after four successive measurements of 9 and ip the success probability has increased to 
approximately 0.93750, while after twelve measurement processes it reaches approximately the value 0.99976. The 
fact that the success probability is maximized for 9 = ir/A indicates that in this case the 9 and ip observables define 
two mutually unbiased bases in a two-dimensional Hilbert space. Fig. QJd) shows p St N versus N for 9 equal to 
7r/12 (circle), 7r/8 (square), and 7r/4 (triangle). Since mutually unbiased bases give the optimal process for each N, 
Ps ,N approaches to 1 faster than in the other cases. As is apparent from Fig. (Dd), the convergence of the success 
probability strongly depends on the relation between the involved bases. In the following section we study this relation 
and generalize the results of this section to the d-dimensional case. 



III. D-DIMENSIONAL CASE 



Two noncommuting and nondegenerate observables defined on a d-dimensional Hilbert space can have at most n 
equal eigenstates with < n < d — 2. Thus, on this n-dimensional subspace both observables can be well defined 
simultaneously, that is, the system can be described by one of the n common eigenstates. However, on the {d — n)- 
dimensional subspace only one of them can be well defined. This property for two noncommuting observables turns 
up in a Hilbert space only when its dimension is higher that 2. One can easily conclude that observables having some 
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FIG. 1: Behavior of: (a) p s , ma x as a function of N for three values of {<p\p\ip): 2/3(triangle), 1/3 (square), (circle), (b) p s ,jv 
as a function of N for three values of 9: 7r/4(triangle), tt/8 (square), n/12 (circle), with (ip±\p\ip±) = 1. 



common eigendirections are not useful for our purpose. In this case the noncommutativity of the observables is a 
necessary but not sufficient condition, as in the two-dimensional case studied in the previous section. This motivates 
us to study the scheme of driving a quantum state by means of measurement in the general case of a d-dimensional 
Hilbert space. 

We now generalize the previous results to the case of a target state, \<px), belonging to a ci-dimensional Hilbert 
space. Let {|<Pi), |<^2)> ■ ■ ■ > \fd)} an d {|1), |2), . . . , \d}} be orthonormal bases defined by the spectral decompositions 
of the non-degenerate observables (p and 9 respectively. Initially the system is described by the p state to be mapped 
onto the known pure state The probability of success after N measurement processes of 9 followed by ip, is given 
by 

d N k d d 

Ps,N = (<pi\p\<pi) + Eton**) E [ II ( E • E 

Pi,jiPh,j2 ■ ■ -PjnA . (8) 

i=2 k=l n=l j"i=2 j n =2 

where pkj is defined as 

i=l 

The process which maps the \<fk) state with k ^ 1 onto the \tpi) state is fundamental in this protocol because it is 
repeated when we do not succeed. The success probability of this process is pk,i- This can be seen as an inner product 
between A\ vectors whose d real, non-negative components are |(i|(,2fc)| 2 (i = 1, ...,<£), that is, pk,i = Y^i=i ^A}. 
This product is maximum when both vectors are parallel, which implies that = ctkA\ for all k. Since these vectors 
are real and the sum of their components is unitary, we deduce that = 1 for all k and that |(i|(^fc)| 2 = | (i\ipk'} | 2 - 
Therefore, we conclude that |(i|<pfe)| = 1/y/d V i, k. This property indicates that, in the optimal case, the two 9 and 
if observables define two mutually unbiased bases. An alternative proof can be obtained by noting that it suffices to 
optimize the first step. That is, we need to project a |<^j_) state onto some element of the {\i)} basis in order to take 
the state out from the subspace orthogonal to the desired direction. The resulting density matrix of this process is 

p=T,\<i\<p±)\ 2 \m (io) 

i 

Now we look for the basis which leads to the state nearest to the target state \ This can be quantified by means 
of the Hilbert-Schmidt distance |l5j . In this case, we need to minimize the expression 

^ = min||p-|^ 1 )(^ 1 ||| 2 . (11) 

Considering the p state, this expression becomes 



d = 2(i-EK%^>I 2 IM>I 2 )- 

i 



(12) 
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Taking into account the properties of the above probabilities, it is clear that the minimum distance is reached under 
the condition 

\{i\v±)\ 2 = Kvi\i)\ 2 , (13) 

which means that the {\i}} basis must be complementary to the original {l^j)} basis. That is, the two required 
bases are related by means of the discrete Fourier transformation. In this scheme only two complementary bases are 
required, which can always be found 0,0,0]. A different proof can be obtained by interpreting pk,i as a correlation 
function and considering the property p\ > pk [Tol ] . 

For mutually unbiased bases the success probability, Eq. , simplifies considerably to 



P.,JV = l-(l-<Pi|p|¥>i>)(l- 5 J • (14) 
In the limit, d> 1, this expression becomes 

P.,w = l-(l-fai|/»l¥>i»e-3. (15) 

Thus, in the case of higher dimensions, in order to reach a success probability close to 1, it is required that the number 
TV of measurement processes of the 9 observable followed by (p must be larger than the dimension d of the Hilbert 
space. Otherwise, the term (ipi\p\tpi) entering in Eqs. i|14|) and i|15[l dominates. 

We now proceed to obtain an average success probability which does not depend on the initial pure state. This is 
achieved by integrating over the whole Hilbert space, that is 

Ps,N = J dipp s ,N, (16) 

where dtp denotes the Haar integration measure and we consider initially pure states only. In this case the average 
probability p s JV is 

p s ,N = i - (i - + (i - b N I #KviMI 2 , (17) 



d' v d 

where we have considered the case of mutually unbiased bases. The starting point is the identity j2Fj| 

where {\n}} with n = 1, . . . , d is an arbitrary base for a d-dimensional Hilbert space. Taking the trace of this identity 
we obtain 

#<V>|n)<A#> = (19) 
Making n — k and considering \tpi) as belonging to the basis, we obtain, for any \ipi) state, the identity 

<fVltei|V>l a = ^ (20) 

Thereby, the average success probability becomes 

P a , N = l-(l-\) N+1 . (21) 

Thus, if we randomly select an initial state, the average success probability of mapping this state onto the target state 
\cpi) is given by p s N . In the limit of large N, p s N becomes 

p^^l-e-^. (22) 
These results are equal to the case when the initial state p is I/d, see Eqs. <|14|) and (|15fl . 
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An interesting application of this result arises when we study the case of a target state \ipi) belonging to a bipartite 
system, each system being described by a d-dimensional Hilbcrt space. In particular 

|^ 1 )=a|V)|V±)+^±)lV')- (23) 
Assuming a factorized initial state of the form 

Pl = p® p, (24) 
the success probability of the process which maps this state onto the \<pi) state is given by 

p Sl N = 1 - [1 - (i>\p\m - + 25R(a/3*)| 7 | 2 )](l - 1)", (25) 

where the 7 coefficient gives account of the initial decoherence process affecting the p state. This coefficient relates 
the diagonal coefficients of p to the non-diagonal coefficients through the relation 

= jt/WpW}W\pW), (26) 
with < I7I < 1. It can be shown that the success probability can be upper bounded as 

p SlN < 1 - [1 - 1(1 + 25R(a/3*)| 7 | 2 )](l - l) w . (27) 

Thereby, when 5R(a/3*) > 0, the maximum probability for fixed N is achieved under the condition | 7 | 2 = 1, that is, 
for a pure initial state. However, if "R(aj3*) < 0, the probability is maximum when | 7 | 2 = 0. This means that, for 
states fulfilling the condition ^H(a(3*) < 0, such as the singlet state, the success probability is higher in the case of 
total initial decoherence than in any other case, corresponding the smallest probability to an initially pure state. 

This scheme can also be connected to the application of quantum erasure. If we fix the target state, and consequently 
the 9 and ip operators, then the sequence of measurements will map any initial state onto that same target state. 
Thereby, the overall effect will correspond to probabilistically erasing the information content of the initial state. The 
success probability of this probabilistic erasure will be given by Eq. J3lJ- 



IV. GENERALIZATION TO m ORTHOGONAL TARGET STATES 



The above results can be generalized to the case of m orthogonal target states belonging to a d-dimensional 
Hilbert space. Here we consider again the bases {|^i), \v>2)i ■ ■ ■ , \ fd)} and {|1), |2), . . . , \d}} of the observables <p and 
respectively. Our aim is to map the initial state p onto any of the target states . . . , !</?,„)}. After measuring 

the ip observable and failing, the state of the system is in one of the {\(p m +i), . . . , \<Pd)} states. The probability of 
mapping the system from any one of these states onto any one of the target states by a consecutive measurement of 
the 9 and (p observables is given by 

d m 

p= E ( 28 ) 

j—m-\-l k—1 

This probability can be written as a sum of (d~m)m scalar products of the A 3 - vectors defined in the previous section, 
that is 

d rn d 

p= E EE44- (29) 

The maximum value of this quantity is achieved when each scalar product involves two parallel vectors, that is 

A? — u^)~A\ Vj = m + 1, . . . , d and k = 1, . . . , m and i = 1, . . . , d. (30) 

Since the sum of the elements of each A\ vector is unity, we obtain a^fc ~ 1 Vj, k. Thus, it holds that A\ = A\, that 
is, all the A\ vectors are equal. This implies that any state \i) has the same projection onto all the states belonging 
to the {l^i), 1^2), • • ■ , Wd)} basis, which is possible only if |(</?j|«)| 2 = 1/d. Therefore, the <p and 9 observables define 
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mutually unbiased bases. Considering these types of bases, which optimize the process, the probability of mapping 
the initial state p onto any of m states |</?i), . . . , \<p m ) after N measurement processes of 9 and (p is given by 



l-^(^| /0 |^ fe )Wl-_j . (31) 



k=l 

As an application, let us suppose that we want to generate I copies of each of the m orthogonal states 
1^2), • ■ • , \i> m )i belonging to a d-dimensional Hilbert space. So, we can assume that the m target states belong 
to a multipartite system composed of I identical systems and that they have the form \ipi) = \ipi)i <S> \tpi)2 ® ■ ■ ■ ® \tpi)i, 
where i = 1, . . . , m. Then the probability of generating any of these "state-copies" after N processes of measurement 
of 9 followed by (p is 

N 



Ps,n = 1 - 1 - X>fc|phk>' (l - ^) , (32) 



k=l 



where we have assumed that the initial state is factorized and that each of the I systems is in the p state. Independently 
of the initial condition, this probability is closer to unity when d l 3> m and Nm 3> d l . On the other hand, randomly 
selecting an initially pure state the average success probability of mapping this state onto one of the m target "state- 
copies" is given by the expression 

/ m\ N + 1 

Ps,n = 1 -{ 1 -j) > (33) 

which in the limit d l ^> m behaves as 

p S;JV = 1 - exp ^ v ^ ' j . (34) 

Thus, in this limit the more the number I of copies, the probability of success decreases or converges more slowly to 
1. 



V. CONCLUSIONS 



We have studied a scheme to map an unknown mixed state of a quantum system onto an arbitrary state belonging 
to a set of known pure quantum states. This scheme is based on a sequence of measurements of two noncommuting 
observables. The target states are eigenstates of one of the two observables, while the other observable maps the 
states out of the subspace orthogonal to the one defined by the target states. The success probability turns out to 
be maximal under the condition that the observables define mutually complementary bases. In other words, both 
required bases are always related by a discrete Fourier transform. We have also shown that these results hold in the 
case of arbitrary but finite dimensions. The scheme consists of applying a measurement of the <p observable followed 
by N measurement processes each one composed of 9 followed by tf>, this is, (<p9) N (p. The target states belong to the 
spectral dcscomposition of the (p observable. The success probability quickly converges to unity when the number iV 
of the sequences of measurement processes is larger than the dimension of the Hilbert space. We have connected these 
results to the generation of quantum copies, quantum deleting and pure entangled states generation. The extension 
of these results to the case of continuous variables is under study. 
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